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ABSTRACT 

The development of a supersonic axial-flow compressor has been the subject of 
a limited amount of research over the past 32 years (Ferri, 1956; Klapproth, 
1961; and Savage et al., 1961). During the middle 1970's a supersonic axial- 
flow compressor was constructed, but it encountered a blade failure before 
reaching its design point (Breugelmans , 1975). Many reached the conclusion 
that the supersonic axial-flow compressor was a very difficult, if not practi- 
cally impossible, design problem. However, recent renewed interest in super- 
sonic and hypersonic flight vehicles have rekindled interest in the supersonic 
axial-flow fan. For example, a research project to design, build, and conduct 
experiments on a single-stage supersonic axial-flow fan is now underway at the 
NASA Lewis Research Center (Schmidt et al., 1987; and Wood et al., 1987). 

Although past experimentation on this type of compressor has been rather 
sparse, some useful analytical technology has been developed. One example is 
in the area of aeroelastic stability. Since the aeroelastic stability of the 
NASA supersonic through-flow fan was a concern, an analytical capability was 
needed to predict the unsteady aerodynamic loading. Consequently, a computer 
code based on Lane's (1957) formulation was developed for the case of super- 
sonic axial flow (Ramsey and Kielb, 1987). This presentation will discuss 
this code and its application to the flutter analysis of the NASA Lewis super- 
sonic through-flow fan. 

The flutter analysis was performed by incorporating this code into an existing 
aeroelastic code and applying it to the NASA blade. The analysis (Kielb and 
Ramsey, 1988) predicted the blades to be unstable at supersonic relative veloc- 
ities. As a consequence, the rotor blades were redesigned by reducing the 
aspect ratio to bring the through-flow fan into the stable operating range. 
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UNSTEADY AERODYNAMIC MODEL 


Lane's (1957) formulation for the unsteady pressure distribution was used to 
calculate the unsteady aerodynamic loads. This formulation considers a cascade 
of two-dimensional flat plates with arbitrary stagger (provided the locus of 
blade leading edges is located ahead of the Mach lines) and arbitrary inter- 
blade phase angle. The upper figure shows the cascade geometry, and the lower 
figure defines the airfoil unsteady pressure distribution. 


CASCADE GEOMETRY 



UNSTEADY PRESSURE 
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COMPUTER CODE VERIFICATION - PITCHING MOTION 

The pressure distribution and lift and moment coefficients due to torsional 
motion were compared with Nagashima and Whitehead's (1977) published results. 
Close agreement can be seen. 


PRESSURE DIFFERENCE 

(M = 2.5, s/c = 1.0, STAGGER ANGLE = 60°, K = 1.0, x„ = 0.5) 
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MOMENT COEFFICIENT 

(M = 1.2, s/c = 1.0, STAGGER ANGLE = 0°, x 0 = 0.5) 
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COMPUTER CODE VERIFICATION - PLUNGING MOTION 

It would have been ideal to compare the lift and moment plunging coefficients 
with those of Nishiyama and Kikuchi (1973). However, it was felt that the pub- 
lished graphs were too small to accurately digitize. Therefore, the plunging 
coefficients obtained from this code were compared to those of Jordan (1953) 
for an isolated airfoil in supersonic axial flow. Close agreement can be seen. 
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REDUCED FREQUENCY BASED ON CHORD, K 
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STRUCTURAL MODEL 

The classical typical section is used to model the structure. Each airfoil is 
assumed to be a two-degree-of-f reedom oscillator supported by bending and tor- 
sional springs. The airfoil is assumed to be rigid in the chordwise direction. 
Coupling between bending and torsional motions is modeled through the offset 
distance between the center of gravity and the "elastic axis." 


I 
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ORIGINAL DESIGN - 73.3-PERCENT SPAN 


The NASA blade is much higher in solidity and lower in stagger angle than typi- 
cal fan stages. However, the airfoil cross section is similar to that of con- 
ventional fan blades. The first mode is primarily bending, and the second mode 
is primarily torsion. The physical properties of the 73.3-percent span loca- 
tion were chosen as being representative and were used in the flutter analysis. 


NUMBER OF BLADES 58 

MASS RATIO 456.2 

RADIUS OF GYRATION (MIDCHORD) 0.431 

STAGGER ANGLE, deg 28 

SOLIDITY 3.215 

AT DESIGN POINT 

AXIAL MACH NUMBER 2.300 

RELATIVE MACH NUMBER 2.606 

REDUCED FREQUENCY (BENDING MODE) 0.376 

REDUCED FREQUENCY (TORSION MODE) 0.663 
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FLUTTER ANALYSIS 


The flutter analysis was performed by incorporating this unsteady aerodynamic 
code into an existing aeroelastic code that solves the stability problem. The 
aeroelastic code was then applied to the NASA through-flow fan blade (Kielb 
and Ramsey, 1988). The analysis predicted that the through-flow fan would be 
torsionally unstable at supersonic relative velocities. As a result, the 
blade aspect ratio was reduced in the final design to bring the rotor into the 
stable operating range. 


TORSIONAL FLUTTER 



RELATIVE MACH NUMBER 
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SUMMARY 


Lane's (1957) analytical formulation, of the unsteady pressure distribution on 
an oscillating two-dimensional flat plate cascade in supersonic axial flow, has 
been developed into a computer code. This unsteady aerodynamic code has shown 
good agreement with other published data. This code has also been incorporated 
into an existing aeroelastic code to analyze the NASA Lewis supersonic through- 
flow fan design. A more sophisticated aerodynamic model that takes into 
account blade camber and/or thickness is being considered as a follow-on to 
this work. 


• LANE’S (1957) FORMULATION HAS BEEN DEVELOPED INTO AN UNSTEADY 
AERODYNAMIC CODE 

• THE UNSTEADY AERODYNAMIC CODE HAS SHOWN GOOD AGREEMENT WITH 
PREVIOUSLY PUBLISHED DATA 

• THE UNSTEADY AERODYNAMIC CODE HAS BEEN INCORPORATED INTO AN AEROELASTIC 
CODE 

• AN UNSTEADY AERODYNAMIC MODEL THAT INCLUDES THICKNESS AND/OR CAMBER 
EFFECTS IS BEING CONSIDERED FOR FUTURE WORK 
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APPENDIX - SYMBOLS 


A 

a 0 

a e 

B 

b 

c 

^L^plung 

(ClHots 

(CM)piung 

(CM)tors 

h 

i 

IM() 

Jo 

K 

K h 

L 

M 

M 

P_ 

P + 

RE() 

s 

t 

U 
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2B yi 

speed of sound 
elastic axis position 

Vm 2 - 1 

semi-chord 

chord 

lift coefficient due to plunging motion 
lift coefficient due to pitching motion 
moment coefficient due to plunging motion 
moment coefficient due to pitching motion 
plunging displacement 
imaginary unit 
imaginary part of () 

Bessel function of the first kind of order 0 

reduced frequency based on chord, coc/U 

bending stiffness 

torsional stiffness 

aerodynamic lift 

Mach number 

aerodynamic moment 

pressure on lower surface of airfoil 
pressure on upper surface of airfoil 
real part of () 
blade spacing 
time 

free-stream velocity 
s treamwise coordinate 
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x/c coordinate of pitching axis with respect to the leading edge 

transverse coordinate 

complex amplitude of incidence 

€ n = 1 if n = 0; € n = 2 if n > 1 

dummy variable of integration 

KM/B 2 = wc/B 2 a 0 

air density at free stream 

interblade phase angle 

complex amplitude of dimensionless pressure difference 
o + 2 kMxi 
angular frequency 
unit step function 
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